Direct Measurement of Turbulent Shear 
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A photon correlation method is introduced for measuring components of the shear rate tensor in a 
turbulent soap film. This new scheme, which is also applicable to three-dimensional flows, is shown 
to give the same results as Laser Doppler velocimetry, but with less statistical noise. The technique 
yields the mean shear rate s, its standard deviation a, and a simple mathematical transform of the 
probability density function P(s) of the shear rate itself. 
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I. INTRODUCTION 

Fluids dissipate energy as they flow through pipes or 
past any smooth or rough surface. Examples include river 
flow or wind blowing across the land. This energy dissi- 
pation is proportional to the velocity gradient, or shear 
rate, of the flow at the bounding surface. This frictional 
energy loss, and its dependence on the Reynolds number 
of the flow [l], Q , is not yet fully understood a century 
after the first explanation was advanced by L. Prandtl 
1. 

Here we introduce a new scheme for measuring the 
shear rate near a bounding surface. It also might be 
applicable in the interior of a fluid. Unlike some widely- 
used methods [l|, the shear rate s is recorded at a single 
"point" of size w. The motivation for developing this 
technique was to improve the usual method for measuring 
the shear rate in turbulent flows @, HJ ■ 

The scheme introduced here is that of photon corre- 
lation spectroscopy (PCS) 5]. It is a variant of that 
used by Fuller and Leal to study laminar flows @. For 
turbulence, the shear rate s is a random variable. The 
PCS method enables determination of the time-averaged 
shear rate s, its standard deviation cr, and the gaussian 
transform of the probability density function (PDF) P(s) 
itself. Because the method has not been used before, the 
values of the mean shear s obtained by PCS are com- 
pared with those measured by laser Doppler velocimetry 
(LDV) 0. 

The PCS scheme has the advantage of improved signal- 
to-noise, short data-collection times, and also the com- 
pactness of the apparatus. The PCS scheme can be used 
when the mean flow rate is absent or present. Hence it 
may be useful outside of the domain of turbulence stud- 
ies. With the PCS scheme, a single beam illuminates a 
group of moving particles that scatter light into a pho- 
todetector at some scattering angle 9. The inset of FigfJJi 
shows the incident and scattered laser beam of momen- 



tum k and k s , respectively and the scattering vector 
k = k s — ko. At a point in the flowing soap film, an in- 
cident beam is focused to a bright spot of size w. The 
intensity Iq of the incident beam is taken to be gaussian 
form, Io( r ) = l(0)e~( r l w \ Figured]}, a side view of the 
setup, will be discussed below. 

The velocity at any point r can then be written as the 
velocity at the center of the spot r—0 plus a term propor- 
tional to the shear rate tensor S, which is the quantity 
of interest. The dominant component of S near a wall in 
this experiment is s = d y u(y,t), where u is in the flow 
direction x while y is in the transverse direction in the 
film plane. Note that s is a scalar quantity. Let u(t) 
be the velocity of an illuminated particle at a horizontal 
distance y from the center of the incident beam (y — 0). 
Then 



u{y,t)=u(Q,t) + sy(t) 



(1) 



where the higher order terms have been neglected. 

Within a multiplicative constant, the scattered electric 
field from N particles within the incident beam at time 
t is 



N N 
3=1 



2s(k-r_j )t 



(2) 



Here -Eo( r j) is the incident gaussian field at the posi- 
tion of the j th particle. Because the scattering from 
micron-size particles is almost perfectly elastic, k — 
(Ann/ A) sin(0/2) where A is the vacuum wavelength of 
the incident light beam (633 nm) and n is the refractive 
index of the soap film, which is 99 % water. 

It will first be assumed that the flow is laminar, so 
that S is time- independent, that is to say, the PDF of 
the shear tensor is a delta function centered at the mean 
value of S. Then the intensity correlation function which 
is simply related to the electric field autocorrelation func- 
tion through the Bloch-Siegert theorem [j| (which is ap- 
plicable to any gaussian PDF, including a delta function) 
is g(r) = (I(t)I(t + T))/(I(t)) 2 = 1 + G(r), where 



G(T) = \(E(t)E*(t + T))\ 2 /(I(t)y< 



(3) 
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Evaluating (3), using (l)-(2) and averaging over t gives 
a result previously obtained by Fuller et al. @ for lami- 
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nar flow, as opposed to a turbulent one. They evaluated 
G(t) rather than <?(t). In the experiments described be- 
low, the turbulent soap film flows in the x direction with 
mean velocity U , where this average is over the width W 
of the soap film. Then G(r) = e~ k w s ' r / 2 , where s t 
is an average over time. Use has been made here of the 
gaussian form of the incident beam. 

Because s is a random function for turbulent flows, an 
additional average over s is needed, giving 



G(r) 



2 „2„,,2 2 



/2 P{s)da 



(4) 



(with P(s) having its maximum near s); G(r) is the gaus- 
sian transform of P(s). 

Two important parameters, in addition to w, are W 
and the Reynolds number, Re — UW/v, where v is the 
kinematic viscosity of the soap solution. 

If the supporting walls that bound the film flow are not 
smooth, their roughness R* is another important control 
parameter. As in three dimensions one expects [il-fToj the 
dimensionless factional drag / = vs/U 2 to be indepen- 
dent of Re when it is sufficiently large. It now depends on 
the ratio R* /W At intermediate values of i?e, experi- 
ment and theory [h| support the result / = Gi?e -1 / 2 , 
where C is just a number, and in 3D flows, / = G-Re -1 / 4 . 

If the shear tensor has more than one component @ 



G(r) 



-2Dk t — U t 



A" 2 / e^^ W ' 2 P( S )d S , (5) 



with Sij = S 



jji when the fluid is incompressible, as in this 
experiment [1J . The factors to the left of the integral take 
into account the extraneous effects of particle diffusion 
and transit time broadening; they are discussed later in 
the text. 




Photo Detector 



FIG. 1: (a): Setup for vertically flowing soap film. The film 
flows down from reservoir RT through valve V between strips 
RW and SW, separated by width W. The weight W keeps 
the the nylon wires taut. Inset shows scattering diagram, 
(b): Side view of the setup, showing laser source, focusing 
lens and photodetector 



II. THE EXPERIMENTS 

The experiments were performed on a soap film chan- 
nel, shown in FigJTJwith W=2 cm. The flow is driven by 
gravity, but there is an appreciable opposing force from 
air friction. However, near the vertical plastic strips that 
support the film, the viscous force from the wires domi- 
nates [2J. 

These strips are glued to thin plastic wires 0.5 mm 
in diameter that join to form an inverted V at the top 
and at the bottom, as indicated in Fig[TJ At the apex, 
a small tube connects the reservoir to a valve V that 
controls the flow rate. The wires at the bottom connect 
and deliver the spent soap solution to reservoir RB, where 
it is pumped back to the top reservoir RT to keep the flow 
rate steady. Typical flow rates of the soap solution are 
- 0.2 ml/s. 

The soap solution is 1% Dawn dishwashing detergent 
in water. It is loaded with neutrally buoyant polystyrene 
particles which scatter the incident beam from a 5 mW 
633 nm He-Ne laser into the photodetector, a Perkin 



Elmer SPCM-AQR-12-FC. The laser source is located 
behind the soap film while the photodetector is located 
in front of it as shown in FigQJi. The laser beam is fo- 
cused onto the soap film with a lens of focal length 25 
cm. The photon stream is delivered to the photodetector 
through an optical fiber, where the fiber tip is located 7 
cm from the illuminated spot on the film. 

In these experiments, w is limited by the wavelength 
of visible light, focal length of the focusing lens and the 
diameter of the incident beam and has a value of w = 
100 /im. The scattering vector k is in the vertical (x) 
direction and the dominant component of the shear rate 
is s = dyuijjjt). The diameter <fi of the seed particles 
(0.4 /im) is sufficiently small that their Stokes number 
in the strongest turbulence is less than 0.1 |2[. Hence 
the particle velocities are adequately close to that of the 
fluid. The refractive index of the soap solution is roughly 
1.3. Typically, the scattering angle 8 = 35°, k = 6 x 10 6 
m -1 . Using a seed-particle density of 1.5 gm/1 yields an 
average photon counting rate of 10 6 Hz. 
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Experiments were performed with a horizontally ori- 
ented comb penetrating the soap film at a point above the 
measuring point and with the comb absent. Only with 
the comb present is the turbulence reasonably developed 
and the ener gy sp ectrum is of scaling form, E(k) oc 
with b ~ 3 (2I . Illlfl3 |. This is the enstrophy range, defined 
as the interval where vorticity of larger size fluctuations 
cascade to smaller scales. In two dimensions there is also 
a cascade of energy fluctuations to larger scales, where b 
= 5/3, as in three dimensions. However, it is not acces- 
sible for decaying turbulence, as in this experiment [T^ |. 
By making the bounding walls rough, so that turbulence 
is constantly being generated there, the inverse energy 
cascade can also be seen The teeth of the comb as 
well as their spacing is 2 mm. 

To further test the PCS technique, measurements are 
also made with the comb absent. In this case, there is no 
well-defined energy spectrum decaying as a power law. 
Nevertheless the flow is far from laminar, so that s can 
be measured by both PCS and LDV. 

To first order, e -k 2 w 2 s 2 r 2 /2 _ 1 _ fc 2 w 2 s 2 r 2/ 2 and 

G(r) — > 1 — k 2 w 2 s 2 T 2 /2. However, experimentally G(r) 
is found to be a non-gaussian function. If P(s) oc 

e ~( S -s) 2 /2a 2 



G(r) 



\Jk 2 w 2 a 2 T 2 + I 



-k 2 w 2 ~^T 2 /(2k 2 w 2 <j 2 T 2 +2) 



> (6) 



which is clearly non-gaussian in r. Both panels of Fig. 
[3] show that while P(s) is close to gaussian form, the 
gaussian fit (solid lines) is not perfect. These "good" fits 
to gaussian form were unexpected. 

There are two other effects that can contribute to the 
decay of G(r): thermal diffusion of the seed particles 
and transit time broadening, which can be dominant for 
large U/w. Both of these contributions are small in these 
experiments but are easy to correct for [l4| • To take dif- 
fusion into account, one multiplies Eq. (4) by the factor, 
Gd = e~ 2Dk r , where D is the diffusion constant, which, 
for spherical particles of diameter <j> is ksT '/3ttt]<P, where 
ks is Boltzmann's constant and 77 is viscosity. 

As for the transit time effect, particles passing through 
a beam of size w produce a burst of light inten- 
sity that temporally modulates the scattered light. 
The multiplicative correction factor here is Gu(t) = 
P -(u 2 ~ 



'AO 



14]. 



The decay times for both of these effects is long com- 
pared to the viscous decay time of interest so these mul- 
tiplicative time factors can be dropped. For example, 
with a spot size w = 100/im and a typical mean velocity 
of U = 1 m/sec, the transit time t u associated with the 
effect is or order a/U ~ 0.1 ms. This is fifty times longer 
than typically measured r c . Diffusion times are much 
longer than this and hence contribute insignificantly to 
the decay of G(r). 

FigH] shows G(r) for measurements made with the 
comb absent (a) and present (b), respectively. Here U ~ 
2 m/s in both experiments, W=2 cm, and w — 100 /im. 



The vertical axis is linear, but the horizontal axis is log r, 
so as to display several decades of lag time. The insets 
to both figures show log G(r) vs r 2 , so that a gaussian 
decay of G(r) appears as a straight line. The straight 
lines in the lower insets indicate that G(r) is indeed of 
gaussian form for very small r. They are a best fit to 
the experimental curves and correspond s = 1600 s _1 
and 1000 s _1 for the experiments with and without the 
comb. 

The solid lines in the upper insets to FigfS] are best 
fits under the assumption of a gaussian P(s). A good 
fit clearly extends beyond the small-r limit and enables 
the determination of the standard deviations a of the 
mean shear as well as s itself. The mean shear s is cal- 
culated from the definition of variance, a 2 = s 2 — s 2 ,s = 

V ' s 2 — a 2 . The results are s = 950 Hz a = 300 Hz with 
the comb absent and s — 1620 Hz a = 500 Hz with the 
comb present. The ratio of a to s is near 20 %. 

The shear measurement is done in the viscosity- 
dominated layer of width S(x), where x is the distance 
from the comb. Ideally the spot size w should be much 
smaller than 5(x). The function u(y) is proportional to 
y within 5(x). Prior experiments have established that 
at x = 20 cm below the comb, where the measurements 
were made, 8 is roughly 200 /zm [2|. Thus the beam size 
w is small enough to correctly measure the viscous shear 
rate. 

The single-point PCS measurements of s are now 
compared with those of LDV, made in the traditional 
way; the vertical velocity component u is measured at 
two nearby horizontally-spaced points in the viscosity- 
dominated interval. 

The LDV measurements were made 2.5 cm below the 
PCS beam spot, which is 80 cm below point P in FigQ] 
and 20 cm below the comb. The LDV measurement point 
is advanced in 50 /im steps starting at y—0. The min- 
imum useful value of y is dictated by the necessity of 
avoiding strong light scattering from the supporting plas- 
tic strip with its edge at y—0. 

The LDV laser source is 514 nm line from a Coherent 
argon-ion laser operating at a power of 500 mW, roughly 
one hundred times that used in the PCS device. The data 
collection time for each measurement of u(y) is roughly 
20 s. Because the correlation time is of the order of mi- 
croseconds, and the counting rate is of the order of MHz 
the function form of G(r) emerges after only a few sec- 
onds of data collection with the correlator. 



III. RESULTS 

The mean shear rate s in the viscous region obtained 
by LDV and PCS agree to within one standard deviation, 
as seen in Table [I] The uncertainties are deduced from 
seven measurements made at the indicated values of U. 
From an individual run, one cannot extract a from the 
LDV data, because noise fluctuations can change even 
the sign of d y u(y, t). 
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(a) Normalized Correlation Function Near Wall (No Comb) (b) Normalized Correlation Function Near Wall (Comb) 




Lag time (ms) Lag time (ms) 

FIG. 2: Panel (a) and (b) are semilog plots of typical correlation functions G(t) obtained with the comb absent and present. 
The conditions of these measurements are described in the text. The first data point in panel (a) should be ignored; it is 
instrumental in origin. The curved solid lines are a best fit to the data using a Gaussian P(s). The straight solid lines are first 
order approximation of G(r) and valid only for k 2 w 2 s 2 T 2 << 1 . However, a Gaussian P(s) fits G(r) reasonably well for more 
than a decade. This fit has no theoretical basis. 



TABLE I: Mean shear rate s as measured by LDV and PCS 
in a narrow range of mean flow speeds (comb inserted). 



Shear Rate (Comb) 



U (m/s) 


LDV s (Hz) 


PCS s (Hz) 


1.4 


1500 


1660 


1.6 


1120 


1030 


1.7 


1760 


1830 


1.8 


1880 


1430 


1.9 


1230 


1200 


2.2 


1760 


1700 


2.2 


1860 


1640 


Average 


1590 


1500 


Std. Dev. 


300 (20%) 


300 (20%) 



a iooo 



PCS • 

LDV ▼ 



0.5 1 
Distance from wall (mm) 



The LDV and PCS measurements span the range 
29000 < Re < 45000 and from 40000 < Re < 57000, 
with and without comb respectively. With the comb in 
place, the Taylor microscale Reynolds number Re\ = 
Urms^/v = 130, where A = u rms / a/ ((du/dx) 2 ) = 1 mm. 

The errors from one run to another are not statisti- 
cal in origin. Rather, the source is variations in the 
flow speed through the valve and the motion of the film 
plane caused by velocity fluctuations of the surrounding 
air which could be only partially suppressed by placing 
the entire apparatus in a tent. 

FigEl shows measurements of s as a function of y in 
units of 50 nm obtained using PCS (circles) and LDV 
(triangles) in the range out to y =1.50 mm with the comb 
present. Here U = 2.16 m/s, W=2 cm and the kinematic 
viscosity of the soap solution is close to that of water (y 
=0.01 cm 2 /s), i?e=45,000. 

The main messages conveyed by this graph are (a) the 
two schemes give roughly the same results for s(y), (b) 
near one of the walls, the LDV measurements are noisier 
(for reasons already discussed), and (c) s decreases with 



FIG. 3: Plot of mean shear rate s as a function of distance 
from the wall (in mm) with a comb in place to strengthen the 
turbulence. The mean flow speed U = 2.16 m/s, J?e=45,000. 



increasing y. Even in the absence of air friction, this 
decrease is expected and is well-studied in 3D flows [l|. 
In soap film flows, air friction slows the flow far from the 
walls, making analysis of the data there difficult. 

This experiment indicates it should be possible to mea- 
sure s near the wall and in the interior of 3D flows, though 
care must be taken to collect scattered photons from only 
a small volume in the fluid. Far from a bounding wall in 
3D turbulence, the PCS method will suffer from the lim- 
itation that w should be smaller than the smallest eddy 
size [lj|, defined as r\ = (^ 3 /e) 1 / 4 . Even in these soap 
film experiments, r\ is estimated to be comparable to or 
smaller than w. Yet, as Fig |3] shows, the LDV measure- 
ments of s agree with the PCS result up to y = 1.5 mm 
from a wall, well outside the viscous region. 
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IV. SUMMARY 

Though the photon correlation scheme has been used 
here to measure properties of the shear rate in a two- 
dimensional soap film, it can be used in three dimen- 
sional flows as well. The PCS method has good signal- 
to-noise, is compact, and uses a laser in the mW range. 
The method yields the variance of the shear rate as well 
as its mean value. The correlation function itself is the 
gaussian transform of the probability density function, 



P(s). 
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